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We generalize the standard quantum adiabatic approximation to tlie case of open quantum sys- 
tems. We define the adiabatic limit of an open quantum system as the regime in which its dynamical 
superoperator can be decomposed in terms of independently evolving Jordan blocks. We then es- 
tablish validity and invalidity conditions for this approximation and discuss their applicability to 
superoperators changing slowly in time. As an example, the adiabatic evolution of a two-level open 
system is analyzed. 
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I. INTRODUCTION 

The adiabatic theorem 0, |3j 0] is one of the oldest and 
most useful general tools in quantum mechanics. The 
theorem posits, roughly, that if a state is an instanta- 
neous eigenstate of a sufficiently slowly varying Hamil- 
tonian H at one time, then it will remain an eigenstate 
at later times, while its eigenenergy evolves continuously. 
Its role in the study of slowly varying quantum mechani- 
cal systems spans a vast array of fields and applications, 
such as energy-level crossings in molecules |4, ^] , quan- 
tum field theory @, and geometric phases recent 
years, geometric phases have been proposed to perform 
quantum information processing llOl lllj , with adia- 
baticity assumed in a number of schemes for geometric 
quantum computation (e.g., 0,0,0,01). Moreover, 
additional interest in adiabatic processes has arisen in 
connection with the concept of adiabatic quantum com- 
puting, in which slowly varying Hamiltonians appear as 
a promising mechanism for the design of new quantum 
algorithms and even as an alternative to the conventional 
quantum circuit model of quantum computation 16, 17j. 

Remarkably, the notion of adiabaticity does not appear 
to have been extended in a systematic manner to the 
arena of open quantum systems, i.e. , qu antum systems 
coupled to an external environment [l8l |. Such systems 
are of fundamental interest, as the notion of a closed sys- 
tem is always an idealization and approximation. This 
issue is particularly important in the context of quan- 
tum information processing, where environment-induced 
decoherence is viewed as a fundamental obstacle on the 
path to the construction of quantum information proces- 
sors (e.g., 01). 

The aim of this work is to systematically generalize 
the concept of adiabatic evolution to the realm of open 
quantum systems. Formally, an open quantum system 
is described as follows. Consider a quantum system 5* 
coupled to an environment, or bath B (with respective 



Hilbert spaces 'Hs,'Hb), evolving unitarily under the to- 
tal system-bath Hamiltonian Hsb- The exact system 
dynamics is given by tracing over the bath degrees of 
freedom [l3 



p{t)=TrB[Uit)psBmHt)l 



(1) 



where p{t) is the system state, pss(O) = p(0) (8) pb(0) is 
the initially uncorrelated system-bath state, and U (t) = 
Texp[— i HsB{t')dt'\ (T denotes time-ordering; we set 
h = 1). Such an evolution is completely positive and 
trace preserving [TsL 0, 0| . Under certain approxima- 
tions, it is possible to convert Eq. (Q) into the convolu- 
tionless form 



p(t) = C{t)p{t) 
An important example is 
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(2) 



1 ^ 

p{t) = -i[H{t),p{t)] + - Y,i[T^it),p{t)T\{t)] 



+ [Ut)p{t),T\{i^] 



(3) 
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Here Hit) is the time-dependent effective Hamiltonian of 
the open system and Fj {t) are time-dependent operators 
describing the system-bath interaction. In the literature, 
Eq. |j2Jl with time-mdependent operators F^ is usually 
referred to as the Markovian dynamical semi grou p, or 
Lindblad equation T?, "27, "25, "23] [see also Ref.^ for a 
simple derivation of Eq. from Eq. (^]. However, the 
case with time-dependent coefficients is also permissible 
under certain restrictions j23] . The Lindblad equation re- 
quires the assumption of a Markovian bath with vanish- 
ing correlation time. Equation Q can be more general; 
for example, it applies to the case of non-Markovian con- 
volutionless master equations studied in Ref. 0| . In this 
work we will consider the class of convolutionless master 
equations Q. In a slight abuse of nomenclature, we will 
henceforth refer to the time-dependent generator C{t) as 
the Lindblad superoperator, and the Ti{t) as Lindblad 
operators. 

Returning to the problem of adiabatic evolution, con- 
ceptually, the difficulty in the transition from closed to 
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open systems is that the notion of Hamiltonian eigen- 
states is lost, since the Lindblad superoperator - the 
generahzation of the Hamiltonian - cannot in general be 
diagonalized. It is then not a priori clear what should 
take the place of the adiabatic eigenstates. Our key in- 
sight in resolving this difficulty is that this role is played 
by adiabatic Jordan blocks of the Lindblad superopera- 
tor. The Jordan canonical form 27], with its associated 
left and right eigenvectors, is in this context the natural 
generalization of the diagonalization of the Hamiltonian. 
Specifically, we show that, for slowly varying Lindblad 
superoperators, the time evolution of the density matrix, 
written in a suitable basis in the state space of linear 
operators, occurs separately in sets of Jordan blocks re- 
lated to each Lindblad eigenvalue. This treatment for 
adiabatic processes in open systems is potentially rather 
attractive as it can simplify the description of the dynam- 
ical problem by breaking down the Lindblad superoper- 
ator into a set of decoupled blocks. In order to clearly 
exemplify this behavior, we analyze a simple two-level 
open system for which the exact solution of the master 
equation (|2J) can be analytically determined. 

The paper is organized as follows. We begin, in Sec. 
im with a review of the standard adiabatic approxima- 
tion for closed quantum systems. In Sec. IIIII we describe 
the general dynamics of open quantum systems, review 
the superoperator formalism, and introduce a strategy to 
find suitable bases in the state space of linear operators. 
Section If VI is devoted to deriving our adiabatic approxi- 
mation, including the conditions for its validity. In Sec. 
IVl we provide a concrete example which illustrates the 
consequences of the adiabatic behavior for systems in the 
presence of decoherence. Finally, we present our conclu- 
sions in Sec. IVII 



II. THE ADIABATIC APPROXIMATION IN 
CLOSED QUANTUM SYSTEMS 

A. Condition on the Hamiltonian 



there corresponds a unique eigenstate, adiabaticity is 
then defined as the regime associated to an independent 
evolution of the instantaneous eigenvectors of H{t). This 
means that instantaneous eigenstates at one time evolve 
continuously to the corresponding eigenstates at later 
times, and that their corresponding eigenenergies do not 
cross. In particular, if the system begins its evolution in 
a particular eigenstate \n{0)), then it will evolve to the 
instantaneous eigenstate \n{t)) at a later time t, without 
any transition to other energy levels. In order to obtain 
a general validity condition for adiabatic behavior, let 
us expand \'ip{t)) in terms of the basis of instantaneous 
eigenvectors of H{t), 

D 
n=l 

with a„(i) being complex functions of time. Substitution 
of Eq. into Eq. ^ yields 

+a„|n)) e-^A' ''^'^^-C*') = 0, (7) 

n 

where use has been made of Eq. © . Multiplying Eq. 
by {k{t)\, we have 



where 



^k{t)=E,,{t)~Ek{t). (9) 



A useful expression for {k\h), for k ^ n, can be found by 
taking the time derivative of Eq. ISJ and multiplying the 
resulting expression by which reads 



{k\n) = — [n ^ k). 



9nk 



(10) 



Therefore, Eq. © can be written as 



To facilitate comparison with our later derivation of 
the adiabatic approximation for open systems, let us be- 
gin by reviewing the adiabatic approximation in closed 
quantum systems, subject to unitary evolution. In this 
case, the evolution is governed by the time-dependent 
Schrodinger equation 



H{t)\m)^^m)), 



(4) 



where H{t) denotes the Hamiltonian and \i-'{t)) is a quan- 
tum state in a I?-dimensional Hilbert space. For sim- 
plicity, we assume that the spectrum of H(t) is entirely 
discrete and nondegenerate. Thus we can define an in- 
stantaneous basis of eigenenergies by 



H{t) \n{t)) = E,,{t) \n{t)), 



(5) 



with the set of eigenvectors \n{t)) chosen to be orthonor- 
mal. In this simplest case, where to each energy level 



^a,(fc|fc)-^a„fflMe~^/o'^*'«-(*'). (11) 



gnk 



Adiabatic evolution is ensured if the coefficients ak(t) 
evolve independently from each other, i.e., if their dy- 
namical equations do not couple. As is apparent from 
Eq. lll|l . this requirement is fulfilled by imposing the 
conditions 



max 

0<t<T 



{k\H\n) 



gnk 



< mm \gnk\ 

0<t<T 



(12) 



which serves as an estimate of the validity of the adia- 
batic approximation, where T is the total evolution time. 
Note that the left-hand side of Eq. (|12|l has dimensions 
of frequency and hence must be compared to the relevant 
physical frequency scale, given by the gap gnk lEH^- For 
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a discussion of the adiabatic regime when there is no gap 
in the energy spectrum see Refs. 29, SOJ. In the case of 
a degenerate spectrum of H{t), Eq. ifTIHl holds only for 
eigenstates \k) and \n) for which En ^ Ek- Taking into 
account this modification in Eq. it is not difhcult 

to see that the adiabatic approximation generalizes to 
the statement that each degenerate eigenspace of H{t)^ 
instead of individual eigenvectors, has independent evo- 
lution, whose validity conditions given by Eq. (fT^ are 
to be considered over eigenvectors with distinct energies. 
Thus, in general one can define adiabatic dynamics of 
closed quantum systems as follows: 

Definition II. 1 A closed quantum system is said to un- 
dergo adiabatic dynamics if its Hilbert space can be de- 
composed into decoupled Schrodinger eigenspaces with 
distinct, time- continuous, and noncrossing instantaneous 
eigenvalues of H{t). 

It is conceptually useful to point out that the relation- 
ship between slowly varying Hamiltonians and adiabatic 
behavior, which explicitly appears from Eq. \i2\ . can also 
be demonstrated directly from a simple manipulation of 
the Schrodinger equation: recall that Hit) can be diago- 
nalized by a unitary similarity tranformation 



Hd{t)^U-\t)H{t)U{t), 



(13) 



where Hd(t) denotes the diagonalized Hamiltonian and 
U{t) is a unitary transformation. Multiplying Eq. by 
U~^{t) and using Eq. l(T^ . we obtain 

Hd\^P)d = i\i^)d-iU-^m. (14) 

where \ip)d = U'-^li/j) is the state of the system in the ba- 
sis of eigenvectors of H{t). Upon considering that H{t) 
changes slowly in time, i.e., dH{t)/dt w 0, we may also 
assume that the unitary transformation U{t) and its in- 
verse U~^{t) are slowly varying operators, yielding 



Hc^{t)m))d^^\m)d■ 



(15) 



Thus, since Hd{t) is diagonal, the system evolves sepa- 
rately in each energy sector, ensuring the validity of the 
adiabatic approximation. In our derivation of the condi- 
tion of adiabatic behavior for open systems below, we will 
make use of this semi-intuitive picture in order to moti- 
vate the decomposition of the dynamics into Lindblad- 
Jordan blocks. 



B. Condition on the total evolution time 

The adiabaticity condition can also be given in terms 
of the total evolution time T. We shall consider for sim- 
plicity a nondegenerate H{t); the generalization to the 
degenerate case is possible. Let us then rewrite Eq. l(TT|) 
as follows ^511: 



gnk 

(16) 



where 7fc(t) denotes the Berry's phase associated to 
the state |fc). 



7fc(0 / dt'{k{t')\k{t')). 
Jo 



(17) 



Now let us define a normalized time s through the vari- 
able transformation 



t^sT, 0<s<l. 



(18) 



Then, by performing the change t ^ s in Eq. (|16|l and 
integrating, we obtain 



ak{s) e 



'hk{s) 



ak 



(o)-E/ 



ds- -— e 

9nk{s') 



where 



Fnk{s) = anis) (fc(s)l^^ln(s)) e^'^'-^^). 
ds 



(19) 



(20) 



However, for an adiabatic evolution as defined above, 
the coefficients a„(s) evolve without any mixing, which 
means that a„(s) « a„(0) e*^"^'*^ Therefore, 

Fnkis) = a„(0) (fc(s)|^^|n(s))e-['^''(^)-^"(^)l. (21) 
as 

In order to arrive at a condition on T, it is useful to 
separate out the fast oscillatory part from Eq. (|19|l . Thus, 
the integrand in Eq. (|19() can be rewritten as 



Fnkis') 

7 — r' 

gnkis') 



-iT ds"g„k{s") 



T 



d_ f Fnkjs') ,T r/ ds"a„,,(s") 

ds' \gUs') 



Substitution of Eq. if^ into Eq. (fT^ results in 

afe(s)e"'^'=(^) = 

'F„fc(0) Fnk{s) ^-iT}^ rfs'ff„fc(s') 



(22) 



flfe 



alkio) glkis) 



+ / ds'e-^^So'^^"9^^(^")±El^ 

ds' glkis') 



(23) 



A condition for the adiabatic regime can be obtained 
from Eq. H23() if the integral in the last line vanishes for 
large T. Let us assume that, as T — > oo, the energy dif- 
ference remains nonvanishing. We further assume that 
d{Fnk{s')/gni;{s')}/ds' is integrable on the interval [0 , s] . 
Then it follows from the Riemann-Lebesgue lemma '32l| 
that the integral in the last line of Eq. H23|) vanishes in 
the limit T ^ oo (due to the fast oscillation of the inte- 
grand) 33] . What is left are therefore only the first two 
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terms in the sum over k oi Eq. (|23|l . Thus, a general 
estimate of the time rate at which the adiabatic regime 
is approached can be expressed by 



F 



where 



F - max |a„(0) (fc(s)|^^|n(s)}|, 
o<s<i ds 

g = min |.g„fc(s)| , 

0<s<l 



(24) 



(25) 



with max and min taken over all k and n. A simplifi- 
cation is obtained if the system starts its evolution in a 
particular eigenstate of H{t). Taking the initial state as 
the eigenvector |rn(0)), with am(0) = 1, adiabatic evolu- 
tion occurs if 



T > 



J" 



(26) 



where 



T = max \{k{s)\ ^ } |m(s))| 
o<s<i ds 

g = min \gmk{s) \ ■ 

0<s<l 



(27) 



Equation (|26|l gives an important validity condition for 
the adiabatic approximation, which has been used, e.g., 
to determine the running_^time required by adiabatic 
quantum algorithms (16i.iil7|. 



III. THE DYNAMICS OF OPEN QUANTUM 
SYSTEMS 

In this section, we prepare the mathematical frame- 
work required to derive an adiabatic approximation for 
open quantum systems. Our starting point is the convo- 
lutionless master equation 0. It proves convenient to 
transform to the superoperator formalism, wherein the 
density matrix is represented by a D^-dimensional "co- 
herence vector" 



\p)) = ( Pi P2 ■■■ pd^ y 



(28) 



and the Lindblad superoperator C becomes a {D^ x D"^)- 
dimensional supermatrix . We use the double bracket 
notation to indicate that we are not working in the stan- 
dard Hilbert space of state vectors. Such a representation 
can be generated, e.g., by introducing a basis of Her- 
mitian, trace-orthogonal, and traceless operators [e.g., 
su(I?)], whence the pi are the expansion coefficients of p 
in this basis [23 , with pi the coefficient of / (the identity 
matrix). In this case, the condition Tip^ < 1 corresponds 
to IIIp))!! < 1, p ^ p^ io Pi — Pi, and positive semidefi- 
niteness of p is expressed in terms of inequalities satisfied 
by certain Casimir invariants [e.g., of su{D)] \34L 1351 is^ . 



A simple and well-known example of this procedure is the 
representation of the density operator of a two-level sys- 
tem (qubit) on the Bloch sphere, via p = (/2 + if • o*)/2, 
where — {(Tx, (Jy, cTz) is the vector of Pauli matrices, 
I2 is the 2x2 identity matrix, and if £ ] 
dimensional coherence vector of norm< 1. 
ally, coherence vectors live in Hilbert-Schmidt space: a 
state space of linear operators endowed with an inner 
product that can be defined, for general vectors u and v, 
as 



is a three- 
More gener- 



{U,V) EE {{u\v)) 



1 



Tr (u'v 



(29) 



where TV is a normalization factor. Adjoint elements {{v\ 
in the dual state space are given by row vectors defined as 
the transpose conjugate of \v)): {{v\ = (i;*, tij, I'Jja)- 
A density matrix can then be expressed as a discrete su- 
perposition of states over a complete basis in this vector 
space, with appropriate constraints on the coefficients so 
that the requirements of Hermiticity, positive semidefi- 
niteness, and unit trace of p are observed. Thus, rep- 
resenting the density operator in general as a coherence 
vector, we can rewrite Eq. (j^J in a superoperator lan- 
guage as 



c{t)\p{t)))^\pm, 



(30) 



where C is now a supermatrix. This master equation gen- 
erates nonunitary evolution, since C{t) is non-Hermitian 
and hence generally nondiagonalizable. However, it is 
always possible to obtain an elegant decomposition in 
terms of a block structure, the Jordan canonical form 
[2^. This can be achieved by the similarity transforma- 
tion 



Cj{t)=S-'{t)£.{t)S{t), 



(31) 



where Cj{t) = diag( Ji, J„i) denotes the Jordan form 
of C{t), with Ja representing a Jordan block related to 
an eigenvector whose corresponding eigenvalue is A^, 



/ Aa 1 

A„ 




V 



Aa 1 

• • • A„ / 



(32) 



The number m of Jordan blocks is given by the num- 
ber of linearly independent eigenstates of C{t), with each 
eigenstate associated to a different block Ja- Since C{t) 
is in general non-Hermitian, we generally do not have a 
basis of eigenstates, whence some care is required in or- 
der to find a basis for describing the density operator. 
A systematic procedure for finding a convenient discrete 
vector basis is to start from the instantaneous right and 
left eigenstates of C{t), which are defined by 



{{Qait)\C{t) 



((Qa(0|Aa(t), 



(33) 
(34) 
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where, in our notation, possible degeneracies correspond 
to Aq, — \p, with a ^ (5. In other words, we reserve a 
different index a for each independent eigenvector since 
each eigenvector is in a distinct Jordan block. It can 
immediately be shown from Eqs. H33|l and (|34() that, for 
A„ ^ A/3, we have {{Qa{t)\Vp{t))) = 0. The left and 
right eigenstates can be easily identified when the Lind- 
blad superoperator is in the Jordan form Cj{t). Denot- 
ing \Va{t))),j = S'~^(0 I'Pait))), i-e., the right eigenstate 
of Cj{t) associated to a Jordan block J^,, then Eq. H33|l 
implies that \T-'a{t)))j is time-independent and, after nor- 
malization, is given by 



The additional left vectors are defined as 



/1\ 





(35) 



where only the vector components associated to the Jor- 
dan block Ja are shown, with all the others vanishing. 
In order to have a complete basis we shall define new 
states, which will be chosen so that they preserve the 
block structure of Cj{t). A suitable set of additional 
vectors is 



^((4°^ I 



(i,o,o, ...,o), 

(O, ...,0,1,0 ), (40) 



which imply the following expression for the left basis 
vector {{£^a\t)\ = j{{£i"^\S-^{t) for C{t): 

{{£i'Ht)\m = {{£i^+'Ht)\ + ((£«(t)| A„(t). (41) 

Here we have used the notation j{{£a^°' = j{{Qa\ 
and j((£'a"°''| = 0. A further property following from the 
definition of the right and left vector bases introduced 
here is 

{{£iHt)\v'-'\m ^ = S^0S^^. (42) 



This orthonormality relationship between corresponding 
left and right states will be very useful in our derivation 
below of the conditions for the validity of the adiabatic 
approximation. 



1 





Vo/ 



\v. 







(36) 

where Ua is the dimension of the Jordan block and 
again all the components outside Ja are zero. This simple 
vector structure allows for the derivation of the expres- 
sion 

Cjit) \Vlj!^))j = + Xait) (37) 



with 



\ra))j and \vi~'^))j 



0. The set 
with j — 0, {ua — 1)| can immediately be 
related to a right vector basis for the original C{t) by 
means of the transformation iVa^t))) — S{t)\Va^))j 
which, applied to Eq. (|37|l . yields 

m pi^' w)) = i-D^r'Hm + ut) iv^^Hm. (38) 

Equation H38|l exhibits an important feature of the set 

||I?^-'''(t)))|, namely, it implies that Jordan blocks are 

invariant under the action of the Lindblad superoperator. 
An analogous procedure can be employed to define the 
left eigenbasis. Denoting by j{{Qa{t)\ = {{Qa{t)\S{t) 
the left eigenstate of Cj{t) associated to a Jordan block 
Ja, Eq. H34|) leads to the normalized left vector 



((Qal 1^ = (O, ...,0,l) 



(39) 



IV. THE ADIABATIC APPROXIMATION IN 
OPEN QUANTUM SYSTEMS 

We are now ready to derive our main result: an adia- 
batic approximation for open quantum systems. We do 
this by observing that the Jordan decomposition of C{t) 
[Eq. (|31|) ] allows for a nice generalization of the standard 
quantum adiabatic approximation. We begin by defining 
the adiabatic dynamics of an open system as a general- 
ization of the definition given above for closed quantum 
systems: 

Definition IV. 1 An open quantum system is said 
to undergo adiabatic dynamics if its Hilbert- Schmidt 
space can be decomposed into decoupled Lindblad- Jordan 
eigenspaces with distinct, time- continuous, and noncross- 
ing instantaneous eigenvalues of C{t). 

This definition is a natural extension for open systems 
of the idea of adiabatic behavior. Indeed, in this case the 
master equation (j^l can be decomposed into sectors with 
different and separately evolving Lindblad-Jordan eigen- 
values, and we show below that the condition for this to 
occur is appropriate "slowness" of the Lindblad superop- 
erator. The splitting into Jordan blocks of the Lindblad 
superoperator is achieved through the choice of a basis 
which preserves the Jordan block structure as, for exam- 
ple, the sets of right | P^^' (t))) } and left \^{{£i'\t)\^ 

vectors introduced in Sec. IIIII Such a basis generalizes 
the notion of Schrodinger eigenvectors. 
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A. Intuitive derivation 

Let us first show how the adiabatic Lindblad- Jordan 
blocks arise from a siniple argument, analogous to the one 
presented for the closed case [Eqs. (fT^ - (fT^ ]. Multiplying 
Eq. (|30|l by the similarity transformation matrix S~^{t), 
we obtain 

Cj\p)).j = \p)).j~S-'\p)), (43) 

where we have used Eq. (pTT|l and defined \p))j = S~^\p)). 
Now suppose that C{t), and consequently S(t) and its 
inverse S~^{t), changes slowly in time so that S^^{t) « 
0. Then, from Eq. (|43|l . the adiabatic dynamics of the 
system reads 

£j{t)\pmj = \pmj- (44) 

Equation 1)441) ensures that, choosing an instantaneous 
basis for the density operator p{t) which preserves the 
Jordan block structure, the evolution of p{t) occurs sep- 
arately in adiabatic blocks associated with distinct eigen- 
values of C{t). Of course, the conditions under which the 
approximation S~^(t) w holds must be carefully clari- 
fied. This is the subject of the next two subsections. 

B. Condition on the Lindblad superoperator 

Let us now derive the validity conditions for open- 
system adiabatic dynamics by analyzing the general time 
evolution of a density operator under the master equa- 
tion (|30|l . To this end, we expand the density matrix for 
an arbitrary time t in the instantaneous right eigenbasis 

{\vfm} as 

/3=1 j=0 

where m is the number of Jordan blocks and is the 
dimension of the block Jp. We emphasize that we are 
assuming that there are no eigenvalue crossings in the 
spectrum of the Lindblad superoperator during the evo- 
lution. Requiring then that the density operator Eq. 
evolves under the master equation (|30|l and making use 
of Eq. (I2HIl, we obtain 

T.T.^f H"^^» + A,|p?))) = 

(3=1 j=l 

E E (46) 

Equation (|46|l multiplied by the left eigenstate ((fa'*! re- 
sults in 

fW=A„rW+r(^+i)-5: ^ 12??)), (47) 

/3=1 J=0 



with ra°' (t) = 0. Note that the sum over (3 mixes dif- 
ferent Jordan blocks. An analogous situation occurred 
in the closed system case, in Eq. . Similarly to what 
was done there, in order to derive an adiabaticity condi- 
tion we must separate this sum into terms related to the 
eigenvalue Aq of C(t) and terms involving mixing with 
eigenvalues Xp ^ Aq. In this latter case, an expression 

can be found for ((£q'-'|P?)) as follows: taking the time 

derivative of Eq. and multiplying by ((^q^I we ob- 
tain, after using Eqs. iglj and 

((4^^|i>?)) = — (((4^) I Zip?)) 

-|-((4^+i)|I>J/')))-((5«|i)J^-^)))), (48) 
where we have defined 

iOpc.it) = Xp{t) ^ Xa{t) (49) 

and assumed A^ ^ A^. Note that, while ajp^ plays a 
role analogous to that of the energy difference gnk in 
the closed case [Eq. ujpa may be complex. A simi- 
lar procedure can generate expressions for all the terms 
{{£i'''\t>'-J~'''')), with k = 0, Thus, an iteration of 

Eq. igHl) yields 

+ ((4^+'^|I>J''"'^)))- (50) 

From a second recursive iteration, now for the term 
{{Si^^^^iv'p^'^^)) in Eq. we obtain 

ERE "'°_„J^.. ^ (^') 

p=l ^9=1 fe,=0 / [ '^J ^pa 

where 

^.-E^^ ' n E =E--- E ' (52) 

s=l ^9=1 ^5=0 j fci=0 fcp=0 

with So = 0. We can now split Eq. 1)47)1 into diagonal 
and off-diagonal terms 

n^-l 

^ E E (53) 

where the terms ^ |I>?)), for Xp ^ Aq, are given by 
Eq. ()51)l . In accordance with our definition of adiabatic- 
ity above, the adiabatic regime is obtained when the sum 
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in the second line is negligible. Summarizing, by intro- 
ducing the normalized time s defined by Eq. I|18|) . we 
thus find the following from Eqs. (|^ and IHSJ- 

Theorem IV. 2 A sufficient condition for open quantum 
system adiabatic dynamics as given in Definition MV. 1\ is: 



max 

0<s<l 



«1, 



E n E 

p=l \q=l /c,=0 



P+Sp 



(54) 



with Xp ^ Aq, and for arbitrary indices i and j associated 
to the Jordan blocks a and (3, respectively. 

The condition (|54ll ensures the absence of mixing 

of coefficients related to distinct eigenvalues Aq in 
Eq. 153|) . which in turn guarantees that sets of Jordan 
blocks belonging to different eigenvalues of £{t) have in- 
dependent evolution. Thus the accuracy of the adiabatic 
approximation can be estimated by the computation of 
the time derivative of the Lindblad superoperator acting 
on right and left vectors. Equation (|54|l can be simplified 
by considering the term with maximum absolute value, 
which results in: 

Corollary IV. 3 A sufficient condition for open quan- 
tum system adiabatic dynamics is 



A/^""" max 

0<s<l 



{{£• 
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'■'fJa 



•C 1,(55) 



where the max is taken for any a ^ (3, and over all pos- 
sible values of i & {0, Ua — 1}, J G {0, n/j — 1}, and 
p, with 



E n E II 

p=l \q=l fc,=0 



(56) 



-i+l+j 
1+J 



1 



{n^-i + l+jy. 



1. 



Observe that the factor A/^"""'' defined in Eq. H5t)|l is just 
the number of terms of the sums in Eq. H54|) . We have 
included a superscript n^, even though there is no explicit 
dependence on Ufj, since j £ {0, ...,np ~ 1}. 

Furthermore, an adiabatic condition for a slowly vary- 
ing Lindblad super-operator can directly be obtained 
from Eq. |MJ, yielding the following. 

Corollary IV. 4 A simple sufficient condition for open 
quantum system adiabatic dynamics is C ~ 0. 

Note that this condition is in a sense too strong, since 
it need not be the case that £ is small in general (i.e., for 
all its matrix elements). Indeed, in Sec. we show via 
an example that adiabaticity may occur due to the exact 
vanishing of relevant matrix elements of C. The general 



condition for this to occur is the presence of a dynamical 
symmetry |37l| . 

Let us end this subsection by mentioning that we can 
also write Eq. H54|l in terms of the time variable t in- 
stead of the normalized time s. In this case, the natural 
generalization of Eq. (|54|l is 



max 

0<t<T 



E n E 



<C min \ujBa\ 
o<t<T ^ 



(57) 



Note that, as in the analogous condition 11 2() in the closed 
case, the left-hand side has dimensions of frequency, and 
hence must be compared to the natural frequency scale 
LUfja- However, unlike the closed systems case, where 
Eq. (|12ll can immediately be derived from the time con- 
dition 124II . we cannot prove here that ujf^a is indeed the 
relevant physical scale. Therefore, Eq. (|57|l should be 
regarded as a heuristic criterion. 



C. Condition on the total evolution time 

As mentioned in Sec. IHl for closed systems the rate 
at which the adiabatic regime is approached can be esti- 
mated in terms of the total time of evolution, as shown 
by Eqs. (|24|l and H26|l . We now provide a generalization 
of this estimate for adiabaticity in open systems. 



1. One- dimensional Jordan blocks 

Let us begin by considering the particular case where 
C{t) has only one-dimensional Jordan blocks and each 
eigenvalue corresponds to a single independent eigenvec- 
tor, i.e., Aq. = A/3 a = /3. Bearing these assumptions 
in mind, Eq. H53|) can be rewritten as 

Vet — AqTq Tq {{£M)-J2rp{{£^\Vp)), (58) 

where the upper indices i,j have been removed since we 
are considering only one-dimensional blocks. Moreover, 
for this special case, we have from Eq. H51|l 



{{£c.\C\Vp)) ^ 



(59) 



In order to eliminate the term AqTo, from Eq. H58|l . we 
redefine the variable ra{t) as 



ro.{t)^Po.{t) e/o^= 
which, applied to Eq. H58|l . yields 



(60) 



Pa = -Pa ((£a|i>a)) " E ((^al^/s)) e"^'" : (61) 
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with 



^I3ait) = / dt' LUf3a{t'). 



(62) 



Equation (|?IT|l is very similar to Eq. (fTT|l for closed sys- 
tems, but the fact that JI^q is in general complex-valued 
leads to some important differences, discussed below. We 
next introduce the scaled time s — t/T and integrate the 
resulting expression. Using Eq. H59|l . we then obtain 

Pa{s) = Pa(0)- / ds'pa(s')$a(s') 



(63) 



' , e 
where $a(s) is defined by 



$o(s) = {{£ais)\-\V^{s))) 



and Vpai-s) by 



Vl3a{s) ^Pl3{s) {{£a{s)\ 



dCjs) 
ds 



(64) 



(65) 



The integrand in the last line of Eq. (|63(l can be rear- 
ranged in a similar way to Eq. for the closed case, 
yielding 



1 

T 



d I Vp^ 
ds 1^2^ 



d Vpa 
ds Ujj^ 



(66) 



Therefore, from Eq. 163|) we have 



Pais) = PaiO)" / ds'pa{s')<S^a{s') 

Jo 



(67) 



the following time conditions, obtained for each Jordan 
block a of C{s), are satisfied: 



T 3> max 

0<s<l 



E 



ds 



ds' UJ^JS') 



(68) 



Equation H68|l simplifies in a number of situations. 

• Adiabaticity is guaranteed whenever V/sa vanishes 
for all a ^ f3. An example of this case will be 
provided in Sec. 

• Adiabaticity is similarly guaranteed whenever 
^/3a(s), which can depend on T through pp, van- 
ishes for all a,P such that Re(ri^Q) > and 
does not grow faster, as a function of T, than 
exp(T| Reri^al) for all a, (3 such that 'Re{Q,pa) < 0. 

• When Ke{il.f3a) — and Im(ri^a) ^ the inte- 
gral in inequality H68|l vanishes in the infinite time 
limit due to the Riemann-Lebesgue lemma [s^l, as 
in the closed case discussed before. In this case, 
again, adiabaticity is guaranteed provided pp{s) 
[and hence Vpais)] does not diverge as a function 
of T in the limit T — > oo. 

• When Re(ri^ct) > 0, the adiabatic regime can still 
be reached for large T provided that Pis^s) contains 
a decaying exponential which compensates for the 
growing exponential due to Re(ri^Q). 

• Even if there is an overall growing exponential in 
inequality H68|l . adiabaticity could take place over 
a finite time interval [0, T*] and, afterwards, dis- 
appear. In this case, which would be an exclu- 
sive feature of open systems, the crossover time T* 
would be determined by an inequality of the type 
T > a-f 6exp(cr), with c > 0. The coefficients a, h 
and c are functions of the system-bath interaction. 
Whether the latter inequality can be solved clearly 
depends on the values of a, 6, c, so that a conclusion 
about adiabaticity in this case is model dependent. 



Thus a condition for adiabaticity in terms of the total 
time of evolution can be given by comparing T to the 
terms involving indices fJ ^ a.. This can be formalized as 
follows. 

Proposition IV. 5 Consider an open quantum system 
whose Lindblad superoperator C{s) has the following 
properties: (a) The Jordan decomposition of C{s) is given 
by one- dimensional blocks. (6) Each eigenvalue of C{s) is 
associated to a unique Jordan block. Then the adiabatic 
dynamics in the interval < s < 1 occurs if and only if 



2. General Jordan blocks 

We show now that the hypotheses (a) and [b] can be 
relaxed, providing a generalization of Proposition IIV.5I 
for the case of multidimensional Jordan blocks and Lind- 
blad eigenvalues associated to more than one indepen- 
dent eigenvector. Let us redefine our general coefficient 
ra\t) as 



(69) 
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which, appHed to Eq. (|53|l . yields 



i|A3=A<, j=0 



- E T.p^ii^^^'^W))^'''''- (70) 

The above equation can be rewritten in terms of the 
scaled time s ~ t/T. The integration of the resulting 
expression then reads 



E E / 



, (71) 



where use has been made of Eq. H51|l , with the sum over 
j and p in the last line denoting 



E-E E n E 

j,p j=a p=l \q=l /c,=0 

The function ^'"^^ (s) is defined by 



(72) 



ds 



(73) 



andyjf)(s)by 



(74) 

The term T ds'p'-a^^\s') in the first line of Eq. ((7T|l . 
which was absent in the case of one-dimensional Jordan 
blocks analyzed above, has no effect on adiabaticity, since 
it does not cause any mixing of Jordan blocks. Therefore, 
the analysis can proceed very similarly to the case of one- 
dimensional blocks. Rewriting the integral in the last line 
of Eq. O, as we have done in Eqs. (|23|) and Ht)7|) . and 
imposing the absence of mixing of the eigenvalues Xp ^ 
Xa, i.e., the negligibility of the last line of Eq. (|71|l . we 
find the following general theorem ensuring the adiabatic 
behavior of an open system. 

Theorem IV. 6 Consider an open quantum system gov- 
erned by a Lindblad superoperator C{s). Then the adia- 
batic dynamics in the interval < s < 1 occurs if and 
only if the following time conditions, obtained for each 
coefficient pa\s) , are satisfied: 



T 3> max 

0<s<l 



E E(-i)'^ 



p+Sp + l 
'■'pa 



(0) 



p+Sp + 1 



ds' 0;?+^'' + ^ 



f3a 



is') 



(75) 



Theorem IIV.6I provides a very general condition for 
adiabaticity in open quantum systems. The comments 
made about simplifying circumstances, in the case of one- 
dimensional blocks above, hold here as well. Moreover, a 
simpler sufficient condition can be derived from Eq. H75|) 
by considering the term with maximum absolute value in 
the sum. This procedure leads to the following corollary: 

Corollary IV. 7 A sufficient time condition for the adi- 
abatic regime of an open quantum system governed by a 
Lindblad superoperator C{t) is 



T > max 

^3 0<s<l 



ds 



(76) 



where max is taken over all possible values of the indices 
Xa / Xp, i, j, and p, with 



A 



f3a 



(n/i-l) (n^-i) I p (j-Sq-i) 

- E E E n E II 

/3|Afj7^A„ j=0 p=l \q=l fc,=0 
{Ua + np - i + 1)\ 



{ria - i + ly.npl 



np 



(77) 



where Ap^ denotes the number of Jordan blocks such that 
Xa ^ Xp. 



D. Physical interpretation of the adiabaticity 
condition 

There are various equivalent ways in which to interpret 
the adiabatic theorem for closed quantum systems 3] . A 
particularly useful interpretation follows from Eq. I|26|) : 
the evolution time must be much longer than the ratio 
of the norm of the time derivative of the Hamiltonian to 
the square of the spectral gap. In other words, either 
the Hamiltonian changes slowly, or the spectral gap is 
large, or both. It is tempting to interpret our results in 
a similar fashion, which wc now do. 

The quantity vj';^^\ by Eq. |(73l, plays the role of the 
time derivative of the Lindblad superoperator. However, 
the appearance of exp[TKeilpa{s)] in Eq. H75|) has no 
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analog in the closed-systems case, because the eigenval- 
ues of the Hamiltonian are real, while in the open-systems 
case the eigenvalues of the Lindblad superoperator may 
have imaginary parts. This implies that adiabaticity is a 
phenomenon which is not guaranteed to happen in open 
systems even for very slowly varying interactions. In- 
deed, from Theorems IIV.5I and IIV.6I possible pictures of 
such system evolutions include the decoupling of Jordan 
blocks only over a finite time interval (disappearing af- 
terwards), or even the case of complete absence of decou- 
pling for any time T, which implies no adiabatic evolution 
whatsoever. 

The quantity ujpa, by Eq. (|^ . clearly plays the role of 
the spectral gap in the open-system case. There are two 
noteworthy differences compared to the closed-system 
case. First, the iopa can be complex. This implies that 
the differences in decay rates, and not just in energies, 
play a role in determining the relevant gap for open sys- 
tems. Second, for multidimensional Jordan blocks, the 
terms ajp^ depend on distinct powers for distinct pairs 
/3, a. Thus certain ujpa (those with the higher exponents) 
will play a more dominant role than others. 

The conditions for adiabaticity are best illustrated fur- 
ther via examples, one of which we provide next. 



V. EXAMPLE: THE ADIABATIC EVOLUTION 
OF AN OPEN QUANTUM TWO-LEVEL SYSTEM 

In order to illustrate the consequences of open quan- 
tum system adiabatic dynamics, let us consider a con- 
crete example that is analytically solvable. Suppose a 
quantum two-level system, with internal Hamiltonian 
H = (w/2) CTz, and described by the master equation 101, 
is subjected to two sources of decoherence: spontaneous 
emission Ti{t) ~ ^{t)a- and bit flips T2{t) = j{t)ax, 
where a- = — ioy is the lowering operator. Writing 
the density operator in the basis {/2, fa;, , CTz}, i.e., as 
p = (/2 -KIT* • CT*)/2, Eq. (1201) resuhs in 



m)) 



OJVx 

_4e2 





ujVy — 2e'^Vx 
-2{^^ + e')vy 
- 2(7^ + 2e^)vz 









Vx 











(78) 

where Vx{t), Vy{t), and Vz{t) are real functions providing 
the coordinates of the quantum state \p{t))) on the Bloch 
sphere. The Lindblad superoperator is then given by 



C{t) 



(79) 

In order to exhibit an example that has a nontrivial Jor- 
dan block structure, we now assume 7^=0; (which can in 
practice be obtained by measuring the relaxation rate 7 
and correspondingly adjusting the system frequency lo). 


















-2e2 


—LO 










-2e2 - 272 





-4e2 








-4e2- 



We then have three different eigenvalues for £(<), 
Ai = 0, 

A2 = — 2e2 — 7^ (twofold degenerate) 

A3 = -46^-272, 

which are associated with the following three indepen- 
dent (unnormalized) right eigenvectors: 



/O 




Vi 

(80) 

with /(7, e) = — 1 — (72/26^). Similarly, for the left eigen- 
vectors, we find 





((^ri = (l//(7,e),0,0,0), 
(0,7',-7',0) , 
-l//(7,e), 0,0,1 



(0), 



(81) 



The Jordan form of C(t) can then be written as 
'Cj(t) 



/ 

-2e2-72 


\ 







-4e^ 



27^ 



(82) 

(observe the two-dimensional middle Jordan block) , with 
the transformation matrix leading to the Jordan form 
being 



//(7,e) \ 



S{t) 



\ 



1 7-^ 
1 
1/ 



(83) 



Note that, in our example, each eigenvalue of C(t) is as- 
sociated to a unique Jordan block, since we do not have 
more than one independent eigenvector for each A^. We 
then expect that the adiabatic regime will be character- 
ized by an evolution which can be decomposed by single 
Jordan blocks. In order to show that this is indeed the 
case, let us construct a right and left basis preserving 
the block structure. To this end, we need to introduce a 
right and a left vector for the Jordan block related to the 
eigenvalue A2. As in Eqs. (|^ and (^OJ, we define the 
additional states as 




= 0,1,0,0 



(84) 



We then obtain, after applying the transforma- 
tions |I?«(t))) = S{t)\vi'^))j and ((£^(01 = 
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J ((£2"-' I S ^{t), the right and left vectors 



I 



V 



, ((4°) I = (0,0,1,0) . (85) 



Expanding the coherence vector in the basis | (i))) | , 
as in Eq. (|45|l . the master equation Ij^Of) yields 



/(7,e)rr^+/(7,e)^r=0 



(0) 



(1) 



= - 2e^ 



„(0) 



7 



'1 



JO) 



^4e2 - 27^) 



' 2 ' 

2^ „(o) 
'3 ' 



(86) 



It is immediately apparent from Eq. H8()|) that the block 
related to the eigenvalue A2 is already decoupled from the 
rest. On the other hand, by virtue of the last equation, 
the blocks associated to Ai and A3 are coupled, imply- 
ing a mixing in the evolution of the coefficients rf^ii) 

and rf^ (t) . The role of the adiabaticity will then be the 
suppression of this coupling. We note that in this simple 
example, the coupling between rj"-* (t) and rg"-* {t) would 
in fact also be eliminated by imposing the probability 
conservation condition Trp — 1. However, in order to 
discuss the effects of the adiabatic regime, let us permit 
a general time evolution of all coefficients (i.e., proba- 
bility "leakage") and analyze the adiabatic constraints. 
The validity condition for adiabatic dynamics, given by 
Eq. yields 



(o)i 



Ai — A3 



72 + 2e2 



«|Ai-A3|. (87) 



We first note that we have here the possibility of an adi- 
abatic evolution even without C{t) « in general (i.e., 
for all its matrix elements). Indeed, solving 7^e/e — 77, 
Eq. (|87ll implies that independent evolution in Jordan 
blocks will occur for e{t) oc 7(i). Since /(7,e) = 
— 1 — (72/2e2) is then constant in time, it follows, from 

Eq. H86|l . that r^\t) is constant in time, which in turn 

ensures the decoupling of r (t) and rg"^ (t) . In this case, 
it is a dynamical symmetry (constancy of the ratio of 
magnitudes of the spontaneous emission and bit-flip pro- 
cesses), rather than the general slowness of C{t), that is 
responsible for the adiabatic behavior. The same conclu- 
sion is also obtained from the adiabatic condition (|54|l . 
Of course, Eq. H87|l is automatically satisfied if C{t) 
is slowly varying in time, which means jit) « and 
e{t) sa 0. Assuming this last case, the following solution 
is found: 



(0), 



4'\t)= r^^^(0)t + rr(0) e 



(0), 



= (-2.^-7^)t 



,(1) 
2 

„(0) 



{t) = 



(0)e 



f)(0)e(--^ 



-27^) t 



It is clear that the evolution is independent in the three 
distinct Jordan blocks, with functions ^ (t) belonging to 
different sectors evolving separately. The only mixing is 



{t) and r2^\t), which are components of the 

(0)/ 



between 

the same block. The decoupling of the coefficients r^'' (t) 

and rg^' (t) in the adiabatic limit is exhibited in Fig. ^ 
Observe that the adiabatic behavior is recovered as the 
dependence of e{t) and 7(t) on t becomes negligible. 




1: Parametric evolution of the coefficients rl^\t) and 
r:,-'{t) for < t < 1. The initial conditions are rf'(O) = 



FIG. 

(0)/ 



JO) 



(0) = 1.0 and the decoherence parameters are taken as 
linear functions of time, i.e., e{t) = eo + at and 7(f) = 70 + at, 
with eo = 1-0 and 70 = 0.5. The master equation is solved 
numerically for a > 0. In the adiabatic regime, corresponding 
to a = 0, the evolution of r^i\t) and rf\t) is decoupled, with 



JO) 



{t) = 1 independently of the value of rg"' {t) 



The original coefficients Vx, Vy, and Vz in the Bloch 
sphere basis {I2 , ax, (Jy, (Tz} can be written as combina- 



tions of the functions ri* •* . Equation 1)881) yields 



vx{0) + {vx{0) - vy{0)h' t) e(-2e^-^')* 
VyiO) + (vxiO) - VyiO)hH) e(-2^'-7=)* 

Vz{0) ^ 



Vxit) 

Vy{t) 

with the initial conditions 



= (-4£^-270t 



1 



/(7,e; 



(89) 



vx{0) = r^°\0) + j-'ri'\0), 

Vy{0) 

Vz{0) 



/(7,e) 



+ 4°^(0), 



(90) 



where now rf'\o) = l/f{-j,e) has been imposed in or- 
der to satisfy the Tvp = 1 normalization condition. The 
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Bloch sphere is then characterized by an asymptotic de- 
cay of the Bloch coordinates and Vy , with approach- 
ing the constant value l/f{-f,e). 

Finally, let us comment on the analysis of adiabatic- 
ity in terms of the conditions derived in Sec. IIV CI for 
the total time of evolution. Looking at the matrix el- 
ements of C(t), it can be shown that, for [3 ^ a, the 
only term V^^^^ defined by Eq. H74|l which can be a pri- 
ori nonvanishing is Vis- Therefore, we have to consider 
the energy difference wia = 4e^ -I- 27^. Assuming that 
the decoherence parameters e and 7 are nonvanishing, 
we have ^13 > and hence ilia > 0. This signals the 
breakdown of adiabaticity, unless V13 — 0. However, as 

we saw above, V13 oc {{£'f''\ t = 27^e/e — 277 and 

thus V13 = indeed implies the adiabaticity condition 
e(i) oc 7(i), in agreement with the results obtained from 
Theorem IIV. 21 In this (dynamical symmetry) case adia- 
baticity holds exactly, while if e(t) is not proportional to 
7(t), then there can be no adiabatic evolution. Thus, the 
present example, despite nicely illustrating our concept 
of adiabaticity in open systems, does not present us with 
the opportunity to derive a nontrivial condition on T; 
such more general examples will be discussed in a future 
publication. 



which are preserved under the adiabatic dynamics. Our 
key results are summarized in Theorems IIV. 21 and IIV.6I 
which state sufficient (and necessary in the case of The- 
orem ^^21 conditions for adiabaticity in open quantum 
systems. In particular. Theorem IIV. 61 also provides the 
condition for breakdown of the adiabatic evolution. This 
feature has no analog in the more restricted case of closed 
quantum systems. It follows here from the fact that the 
Jordan eigenvalues of the dynamical superoperator - the 
generalization of the real eigenvalues of a Hamiltonian 
- can have an imaginary part, which can lead to un- 
avoidable transitions between Jordan blocks. It is worth 
mentioning that all of our results have been derived con- 
sidering systems exhibiting gaps in the Lindblad eigen- 
value spectrum. It would be interesting to understand 
the notion of adiabaticity when no gaps are available, as 
similarly done for the closed case in Refs. p9ll30j|. More- 
over, two particularly intriguing applications of the the- 
ory presented here are to the study of geometric phases 
in open systems and to quantum adiabatic algorithms, 
both of which have received considerable recent atten- 
tion 0,0,|1EIS3.E3- We leave these as open problems 
for future research. 
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